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We study a chain of magnetic moments exchange coupled to a conventional three dimensional 
superconductor. In the normal state the chain orders into a collinear configuration, while in the 
superconducting phase we find that ferromagnetism is unstable to the formation of a magnetic spiral 
state. Beyond weak exchange coupling the spiral wavevector greatly exceeds the inverse supercon¬ 
ducting coherence length as a result of the strong spin-spin interaction mediated through the subgap 
band of Yu-Shiba-Rusinov states. Moreover, the simple spin-spin exchange description breaks down 
as the subgap band crosses the Fermi energy, wherein the spiral phase becomes stabilized by the 
spontaneous opening of a p—wave superconducting gap within the band. This leads to the possibil¬ 
ity of electron-driven topological superconductivity with Majorana boundary modes using magnetic 
atoms on superconducting surfaces. 

PACS numbers: 75.30.Hx, 75.75.-c, 74.20.Mn, 03.67.Lx 


The prospect of performing topological quantum com¬ 
putation [1] has stimulated intense investigations into 
condensed-matter systems harboring Majorana bound 
states [2-13]. One potential platform involves mag¬ 
netic atoms arranged in a regular lattice on an 5—wave 
superconducting substrate [14]. This system has re¬ 
ceived a renewed interest due to recent scanning tunnel¬ 
ing microscopy (STM) data possibly supporting the exis¬ 
tence of Majorana bound states in a self-assembled one¬ 
dimensional (Id) array of atomically-spaced Fe atoms on 
the surface of superconducting Pb [15-17]. 

We focus on the case of an STM-assembled magnetic 
atom chain whose spacing is several substrate lattice sites 
[28], where the overlap of atomic wavefunctions is negli¬ 
gible. In that case each atom acts as an isolated mag¬ 
netic moment giving rise to localized, sub-gap Yu-Shiba- 
Rusinov (YSR) states [20-25] within the superconduc¬ 
tor. Overlap between YSR states leads to the forma¬ 
tion of an effectively spinless sub-gap band that can 
undergo a topological superconducting transition con¬ 
trolled by the magnetic ordering of the magnetic atoms 
(spin chain) [26, 27], see Fig. 1. On the other hand, 
the ordering of the spin chain is dictated by the indi¬ 
rect exchange mediated by electrons in the superconduc¬ 
tor, allowing the possibility of an electron-driven “self- 
organized” topological superconducting phase. 

Although many important aspects of spin chains on su¬ 
perconductors have been theoretically studied by several 
authors [19, 26, 27, 30-37] (see also [38-41]), a generic 
conceptual account of the collective ordering mechanism 
of the magnetic and electronic degrees of freedom in a 
bulk (d > 1) superconductor remains elusive (the Id case 
is special due to 2kF nesting and was studied in [30-32], 
see also [38]). 

In this paper, we establish the existence of spiral mag¬ 
netic order and self-organized topological superconduc¬ 


tivity in the case of a magnetic adatom chain on a 3d 
substrate using an exactly solvable minimal model. We 
provide the conceptual framework to describe the order¬ 
ing mechanism in a general setting, which we find to be 
entirely distinct from the 2kF mechanism that can occur 
for a Id substrate [30-32, 38]. In particular, we show 
that ferromagnetism of the spin chain is unstable to the 
formation of a spiral state due to the presence of sin¬ 
glet superconductivity in the host. The corresponding 
spiral wavevector exhibits a nonmonotonic dependence 
on the spin-lattice exchange coupling to the supercon¬ 
ductor. The wavevector grows rapidly as the exchange 
increases from zero due to a strong indirect spin inter¬ 
action mediated by YSR states, and decreases after the 
YSR band crosses the Fermi energy due to the double¬ 
exchange mechanism. In the latter regime the YSR chain 
exhibits self-sustained p—wave superconductivity with 
the possibility of Majorana boundary modes, depending 
on the lattice spacing. Our main results are summarized 
in Fig. 2, with details related to various technical calcu- 



FIG. 1. (Color online) Magnetic atoms on the surface of a su¬ 
perconductor induce local spin-polarized, subgap YSR states 
(black curves). YSR wavefunction overlap (shaded red) leads 
to hopping and p—wave pairing due to hybridization with sin¬ 
glet Cooper pairs in the host. This hybridization can drive 
the spin chain into a spiral state that, in turn, induces spon¬ 
taneous topological superconductivity within the YSR chain. 
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FIG. 2. (Color online) (a) Exact groundstate phase diagram of Eq. ( 1 ) as a function of lattice spacing and YSR energy e 
(Eq. 5) for integer n (n = 15 shown) and = 2000. Outside the range n+1/4 < kpalT^ ^ n + 3/4 the spin chain is 
an antiferromagnet AE with q = 7r/a, while inside this range it is a spiral (S) whose wavevector is shown schematically in 
gray-scale and plotted in panel (b) along the line kpa = IS.Stt. The topologically nontrivial superconducting (TS^) regions 
are bounded by the red curves, across which either a first order transition to the AE state or a continuous transition into 
the trivial S state occurs. The latter transition occurs approximately along the dotted lines: /A = —{kpa mod =F 7T)lkFCL- 

The system is gapless along the ferromagnetic line e/A — (7r/2 — (kpa modyr)) /kpOL separating the topologically distinct TS^ 
regions with topological invariants ±1. (b) Ground-state wavevector q of the magnetic spiral plotted as a function of £. The 
solid line corresponds to the exact solution of Eq. (1), the long-dashed line to the T—matrix approximation of Eq. (6) and the 
dotted line to the YSR approximation of Eq. (7). The inset shows the magnetic order near the topological superconducting 
(TS^) phases bounded by the vertical dash-dotted lines. 

lations delegated to the Supplemental Material [42]. 

A system of magnetic atoms exchange coupled to a 
3d superconductor may be described by the following 
Bogoliubov-de Gennes Hamiltonian 


order in to arrive at the effective spin Hamiltonian 
= ( 2 ) 


H = + + j^^lSr-CT^r, ( 1 ) 

where ^(k) = ^ ^ electron mo¬ 

mentum, position and mass, respectively, while kp is 
the Fermi momentum), ^ 7 /^ 4 ,,T/^j, —T/^j)^ is the 

4-component Nambu spinor written in terms of electron 
annihilation (creation) operators with spin pro¬ 

jection (T. Here cji and are, respectively, Pauli ma¬ 
trices acting in the spin and particle-hole spaces while 
Sr = S(r — rj)Sj is the sum of spins j at sites Vj 
with lattice spacing a and J is the exchange interaction 
constant. 

We provide the exact phase diagram of Eq. ( 1 ) within 
the classical spin approximation [ 21 ] assuming a planar 
spiral ansatz. Si • Sj = S^co^Qa{i — j), where the spi¬ 
ral wavevector Q is treated as a variational parameter 
determined by minimizing the total electronic energy of 
the system, Qmin = Q- The main results of this proce¬ 
dure are shown in Fig. 2, with technical details provided 
in the Supplemental Material [42]. Below, we focus on 
extracting only the most important and physically rele¬ 
vant ordering mechanisms, which we verify against the 
exact results. 

Magnetic Atom Order — The spin-chain magnetic or¬ 
der is easiest to determine in the limit of small exchange 
coupling, where one may integrate out electrons to second 


The electron-induced spin-spin exchange coupling in 
Eq. (2) is {h = 1) [43-46] 


I{r) oc 


-tArCos{2kFr) + ^ sin^(fcirr) 
zirr'^ r^ 


(3) 


where Ui? is the Fermi velocity and ^ = Ui?/A is the 
coherence length of the superconductor. We note that 
solving for the superconducting order parameter self- 
consistently is not important in the dilute atom regime 
{kp ^ since the YSR energy and electron 

propagator between spins are hardly modified [47]. We 
do not consider the dense atom limit kpa < 1 , where 
a collective depletion of A(r) near the chain should be 
taken into account. 

The first term in square brackets in Eq. (3) is the 
familiar Rudermann-Kittel-Kasuya-Yosida (RKKY) in¬ 
teraction [48-50], while the second, strictly antiferro¬ 
magnetic, term arises from the anomalous component 
of the electronic spin susceptibility and disfavors the 
pair-breaking effect of a polarized exchange field. The 
spin-chain groundstate energy is given by the minimum 
Fourier component of the exchange interaction, Eq = Iq. 

In the normal state (A = 0), one finds from Eq. (3) 
a ferromagnet (g = 0) in the range n + 1/4 < kpajn < 
n + 3/4 with integer n and an antiferromagnet {q = Tija) 
otherwise. We emphasize that the q = 2kp spiral 
phase that may arise when the electron gas is also one¬ 
dimensional [30-32, 38] does not occur here; the spin 
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chain is always in a collinear state. However, supercon¬ 
ducting correlations in a 3d electron gas generally make 
ferromagnetism unstable towards spiral formation due to 
the long-range superconducting correction in Eq. (3) (in 
a Id electron gas the q = 2kp order is hardly modified by 
superconductivity [30-32]). Indeed, for A 7 ^ 0 the mini¬ 
mum Fourier component, ^ — |g|/^, is shifted to the 

wavevector q ^ This magnetic instability is similar 
to the Anderson-Suhl transition in 2 d and 3d spin lat¬ 
tices [44, 51] and results from the compromise between 
the shorter-range ferromagnetic RKKY interaction and 
the longer-range antiferromagnetic interaction. 

The spiral formation is stable and strongly reinforced 
by the inclusion of higher order terms beyond the Born 
approximation in Eq. (3). The crucial effect responsible 
for this behavior is the formation and overlap of the sub¬ 
gap YSR bound states induced by the magnetic atoms. 
Formally, these subgap states are encoded in the single 
magnetic atom T—matrices 


Tj{uj) = JSj ■ a 



° -Y , ( 4 ) 

2 VA2 - w2 ^ / 


where up is the normal state density of states at the 
Fermi energy. The presence of the bound state is signified 
by the subgap pole of det Tj{uj) on the real energy axis 
UJ = ±s, 


^ 1 - {'kJSvfP)'^ 

^ ^ 1 + {ttJSvf/2Y ■ 


( 5 ) 


To capture the modification of Eq. (3) due to the presence 
of the YSR bound states, one may replace the bare mag¬ 
netic atom potential with its T—matrix, Eq. (4). Com¬ 
puting the exchange coupling to second order in T leads 
to [42] 


I{r) = 


e^\ 

A2 ) 2{kFrY 
A‘^cos‘^{kpr) \s 


V p 

- cos(2kpr) 

2pr 


(6) 


2kl 


+ ^ [1 — 3cos{2kpr)] 


where J is expressed through 5 using Eq. (5). While in 
the limit of small J, i.e. 5 ^ A, we recover the per¬ 
turbative result (3), we find for small 5 an enhanced an- 
tiferromagnetic component ^ l/j^j due to the exchange 
coupling mediated via the YSR states [46]. This effect 
strongly increases the spiral wavevector of the spin chain. 
Indeed, in the range n + 1/4 < kpajir < n + 3/4 for inte¬ 
ger n, one finds Iq ^ q^ — \q\^~^A/\s\, with the ground- 
state wavevector q ^ rapidly growing as 5 ap¬ 

proaches zero. This behavior is consistent with the exact 
numerical result, shown in Fig. 2b, as long as \s\/A is 
not too small. The apparent divergence of Eq. (6) as 
l^l ^0 signifies its inapplicability for small enough 5 , 
since the effective exchange energy cannot physically ex¬ 
ceed the splitting between YSR energies caused by their 
finite overlap. 


This suggests that below a characteristic value of 5 , 
a crossover (for two spins) to a new regime occurs. To 
see this we follow [46] and consider two YSR states hy¬ 
bridizing with a Cooper pair from the superconducting 
condensate. This process can be described by the matrix 
element U{r) = Ae~^^^cos{kpr)/{kpr) representing the 
overlap between YSR states, see Fig. 1. The energy dif¬ 
ference between the ground state and the excited state 
with YSR states occupied/unoccupied is 2|5|. As long as 
U < 2\e\, the change of the groundstate energy can be 
computed perturbatively, U^/2\s\^ leading to the second 
term in the brackets of Eq. (6). Thus, for the spin chain, 
Eq. (6) is valid down to the scale jej ^ U{a) ^ A/kpa 
(assuming ^ ^ a), at which point q ^ kpa/^. When the 
Fermi energy crosses the YSR band jej ^U{a) ^ A/kpa, 
the perturbation theory in U/\s\ fails since Cooper pairs 
in the bulk strongly hybridize with the YSR states. This 
signals a fundamental change of the groundstate, and 
corresponds to a p—wave superconducting phase transi¬ 
tion in the YSR chain which harbors Majorana boundary 
states (provided spiral order remains stable). 


We now wish to show that in the regime of strong hy¬ 
bridization (where the effective spin model Eqs. ( 2 ), ( 6 ) 
is inapplicable) ferromagnetism remains unstable to spi¬ 
ral formation (i.e. the topological superconductivity is 
stable). The complication is that a new competing ten¬ 
dency arises once the Fermi energy lies in the YSR band: 
the gain of kinetic energy associated with electrons hop¬ 
ping along the YSR chain is maximized for a ferromag¬ 
netic arrangement of the spin chain since this maximizes 
the YSR bandwidth, analogous to the double-exchange 
mechanism. However, the presence of a spiral leads to 
a superconducting gap in the YSR band, leading to a 
concomitant gain of the corresponding condensation en¬ 
ergy, depicted in Fig. 3. At small g, it is easy to see 
that the energy gain from the opening of the supercon¬ 
ducting gap (^Egap oc — g^(9gA/c^)^ln^ is always larger 
than the loss of kinetic (or RKKY) energy 6E]^i^ oc g^, 
where the coefficient of the logarithm is proportional to 
the square of the YSR gap at the YSR Fermi points ±k^ 
(see Fig. 3) expanded for small g, A/^^ « qdqAj^^ (A/. = 0 
for g = 0 since Cooper pairs cannot hybridize with a fer¬ 
romagnetic chain). Qualitatively, we thus find that fer¬ 
romagnetism is generically unstable to the formation of 
spiral phase with wavevector g oc where the pa¬ 

rameter A oc {dqAk^)‘^ entering the exponent depends on 
5 , kpa, and is provided below Eq. (8). The case g = 0 
can only occur if the gap at the YSR Fermi points, and 
thus A, vanishes for k^ 7 ^ O^n/a. 

We study this possibility and the quantitative variation 
of g on 5 by evaluating the YSR contribution to the total 
energy numerically (the above gap RKKY contribution, 
first term in the brackets of Eq. (6), only weakly depends 









on 1^1 <C A here) 
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= Ek = \l{hu - ef + (7) 

k 

Here the hopping hk = h-k and p—wave pairing A/^ = 
—A_/c energies are the same as those derived in Ref. [26], 
up to the energy shift 5 which enters Eq. (7) only as 
an effective chemical potential for the YSR band. The 
precise forms of A/, are provided in the Supplemental 
Material [42] , but are not essential in the following; their 
salient features are sketched in Fig. 3. In particular, one 
sees that A/^ indeed may have extra nodes away from 
k = due to the long-range nature of the pairing 

(i.e. higher odd harmonics beyond sinka are relevant). 

Generally we find that when = hk^ = 0 have 
simultaneous zeroes double-exchange wins and the sys¬ 
tem exhibits a gapless electronic spectrum and ferromag¬ 
netism in the spin chain. Remarkably, the exact loca¬ 
tion of these gapless critical lines in the phase diagram 
of Fig. 2a can be extracted: 




A 

kpa 


[7r/2 — {kpCi modTr)] . 


(8) 


In the spiral regions of Fig. 2, the magnetic order in the 
close vicinity of Eq. (8) is exponentially suppressed qa ^ 

where A ^ ^ A^kpa ) * ^way from this line q grows 

to the scale q kpa/^ dit e A/kpct^ beyond which it 
crosses over to the T—matrix dependence q^ r\j A/|e|, as 
shown in Fig. 2b for the case kpa = IS.Stt {cc = 0). This 
leads to the non-monotonic behavior of q on the YSR 
energy 5, with a peak around 5 ^ A/kpct of order g^ax ^ 
kpa/^. The magnetic order determined from minimizing 
the electronic energy of Eq. (1), shown in Fig. 2b, indeed 
exhibits this crossover behavior and is well-described by 
the T—matrix approximation when > A/kpa and the 
YSR approximation when kl < A/kpa. 

We now briefly discuss the topological properties of 
the model in Eq. (1), delegating a more detailed ac¬ 
count to the Supplemental Material. As discussed in 
Refs. [34, 35, 54-57], the topological class of the YSR 
chain in the planar spiral configuration is BDI due to hid¬ 
den time-reversal and chiral symmetries. Different topo¬ 
logical phases are therefore distinguished by a Z topolog¬ 
ical invariant (winding number), which can only change 
when the gap in the YSR spectrum closes. This occurs 
when the YSR band crosses the Fermi energy or when 
the system is tuned across the gapless ferromagnetic line 
Eq. (8), leading to the topologically non-trivial and dis¬ 
tinct TS^ phases with winding numbers ±1 (for a fixed 
rotation axis of the spiral), shown in Fig. 2a. As a result, 
the boundary between TS+ and TS“ systems will host 
a pair of zero-energy Majorana bound states protected 
by chiral symmetry, whereas the boundary of either sys¬ 
tem with vacuum will only host a single Majorana bound 



FIG. 3. Spectrum Ek of the YSR band, Eq. (7) For g / 0 
a spectral gap opens and the system gains the condensation 
energy, corresponding to the gray area between the full and 
dashed lines for momenta near i/c*. This occurs only if the 
pairing at these points is non-zero: Ak^ ^ 0. When the 
zeroes of hk and Ak (depicted by dots) coincide, the energy 
(7) is minimized at g = 0, corresponding to a gapless electron 
system and a ferromagnetic spin chain. This phase occurs 
along lines in the phase diagram given by Eq. (8). 

state. The presence of two Majorana bound states also 
occurs for a single system in either of the TS^ phases if 
there is a g ^ —g domain wall [56, 57]. 

Experimental Outlook — The YSR topological gap. At, 
obtains a broad maximum as a function of magnetic atom 
spacing (i.e. from one period n of the topologically non¬ 
trivial part of the phase diagram to the next, kpa 
kpa + tt). This is because if the magnetic atoms are 
spaced too closely (a ^ 0) At oc q^s^^aA/kpa oc Aa/^ 
becomes small, while if the atoms are too dilute (a ^ oo) 
the overall energy scale A/kpa is diminished. The maxi¬ 
mum thus occurs for an optimal spacing Uopt determined 
by gmaxffl ~ 1, giving kpa^-pt ~ At this spacing, for 

kp^ = 2000, we find At ^ 200mK, assuming A = 10 K 
(by comparison for Fe on Pb experiments show At < 1 K 
[15, 16]). 

The influence of anisotropy and temperature on the 
topological phase are expected to be unimportant as long 
as the associated energy scale remains below At. The 
magnetic order can also be destroyed by thermal fluc¬ 
tuations, but such effects will be quenched either by fi¬ 
nite chain size or anisotropy [30-32, 35]. For the case 
of single-ion anisotropy, described by the Hamiltonian 
Hd = we find that when D < 0 (easy-plane) the 

spiral is merely fixed to lie in the x^-plane (i.e. the spiral 
axis lies along z), while for D > 0 the spiral axis lies in 
the plane and is orthogonal to the ^ axis. For suffi¬ 
ciently strong anisotropy D > At ^ 20/ieV we expect the 
groundstate to become a topologically trivial, collinear 
magnetic phase. This restricts the self-organized topo¬ 
logical phase to easy-plane systems or those with suffi¬ 
ciently weak easy-axis anisotropy. 
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Self-Organized Topological Superconductivity in a Yu-Shiba-Rusinov Chain 

Supplemental Material 

In this Supplemental Material we compute the exact groundstate energy of the BdG Hamiltonian as a function of 
the spiral wavevector in Sec. and derive the formula determining the exact subgap YSR spectrum in Sec. , which 
allows us to classify and determine the exact topological phase boundaries in Sec. . 


GROUNDSTATE ENERGY 

We start from the BdG Hamiltonian (Eq. (1) of the main text) and compute the thermodynamic potential Q = 
— ^InZ at inverse temperature [3 using the quantum partition function Z = Tre“^^. This allows us to extrapolate 
to the groundstate energy in the limit /3 ^ oo. Following Ref. [SI] (see p. 254), we introduce the real parameter A 
multiplied onto the exchange part of the BdG Hamiltonian (1), 

^ E ^p + E (SI) 

p j 

written with 4x4 Nambu matrices Jj = Jex{Sj • cr)ro. From this definition, the change of thermodynamic potential 
due to the spin lattice can be found as an integral of dfl/dX over A: 

Afi = fi(l) - m = dX ^ ■ (^2) 

Here we use the matrix notation lij = Jidij and = Q^{ikn] r^), with the Nambu matrix structure being implied. 
The trace involves a sum over spin lattice positions, {i}, Nambu indices and Matsubara frequencies. Notice that the 
Green function is defined with averaging, (.. .)a, with respect to H\. 

In terms of the bare Green function, G^, the dressed Green function matrix is obtained by summing the Dyson 
series 


G^ = GV G^AJG^ + G^AJJG^AJJG^ + ..., (S3) 

We note that the 4-spinor notation, 4^ = exhibits the usual particle-hole redundancy since 4^’t'(x) = 

r^cr^4^(x). The energy overcounting is avoided by the factor of 1/2 multiplying the Hamiltonian (SI). This factor 
of 1/2 also compensates a factor of 2 appearing at each vertex in the Dyson series (S3) due to the fact that the 
perturbation series for the Green function involves both anomalous (T7-4^(r)4^(0))o as well as normal (T7-4^(r)4^’*'(0))o 
Nambu contractions. Since, however, the two are related by multiplication with from the right, their contributions 
are equivalent and give rise to a compensating factor of 2 at each vertex As a result, one arrives at the Dyson 
series (S3) containing only Nambu normal Green functions. 

Inserting (S3) into (S2), the A-integration can be carried out order by order replacing A’^”^ by 1/n, and the infinite 
series can be recognized as a logarithm: 


AD =--Trln 

2/3 


1 - 


(S4) 


Since we are interested in the magnetic ordering arising from the induced interaction between different spins, we 
focus on the part of the energy shift which depends on the spin lattice spacing, a. The only a-independent terms in 
the Dyson equation are those which only involve repeated scattering from a single magnetic atom. To facilitate this 
separation we split the bare Green function matrix into a sum of local (oc 5ij)^ and non-local (oc 1 — 5ij) terms: 


G^ 



(S5) 


The a dependent part of (S4) is therefore found by subtracting the local part, i.e. 


Ai^a = -^(TVln [1-JG°] -IVln [l - JG°J) = 


1 -TG° 


(S6) 
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where the local T-matrix has been introduced as 

T=(1-JGL)”'J. (S7) 

To investigate the magnetic structure, we use a variational approach in which we minimize the electronic ground 
state energy with respect to the magnetic ordering profile. We consider a planar spiral order ansatz whose spiral axis 
of rotation is defined along z 


Sr = S S{x — ja)S{y)S{z) (cos qaj, sin qaj, 0), 

3 


(S8) 


which interpolates between ferromagnetic (F) (g = 0) and anti-ferromagnetic (AF) (g = ±7r/a) order [S2]. Within 
the spiral ansatz, the j dependence of Tj can be conveniently written as Tj = where T = 

J(7x{^ — Since Qq oc ao, the Green function commutes with and we may rewrite Eq. (S6) 

as 


= -^TrlnM, 


.ij = % - (1 - %)go (iK, a{i - j)) 


The local Green function and the T-matrix are given by 

n fl n\ nPFiknTo + AT^ 

rrfj \ T b , iknTo + 


-1 


(S9) 


(510) 

(511) 


As noted in the main text the determinant of T exhibits a pole on the real energy axis ikn e where 5 is defined 
in Eq. (5). Eor r ^ vfI^d^ with Debye frequency ujd representing the high-frequency cut-off in the BCS Green 
function, the non-local Green function takes the form 


Qo{ikn,r) 


FVF e-''V'=^+^V«F 
2 kFV 


iknTp + ATx 


sin(/ci?r)-h T^cos(/ci?r) , 


(S12) 


which applies to propagation along the spin chain, provided kpa ^ 

The matrix M in Eq. (S9) exhibits translational invariance, allowing it to be diagonalized using a set of Bloch 
states labeled by quasi-momentum \k\ < ir/a. In this basis, the operator trace over spin lattice positions in Eq. (S6) 
becomes an integral over the Brillouin zone. Eocusing on the low-temperature limit l3~^ A/Zc^a, we replace the 
Matsubara summation by a frequency integral to obtain the variational groundstate energy for a chain of N spins 


E{q) = 


Na 

~Y 


/ 


^i^Trln 

2f 2f 


1 - Qo{iuj,k- 


(S13) 


In Eq. (S13) we introduced the Eourier transform of the non-local part of the Green function 


gp{iuj,k) 


(S14) 


The remaining trace over the 4x4 spin/particle-hole space in Eq. (S13) can be performed explicitly. The groundstate 
magnetic order is characterized by the wavevector g which minimizes Eq. (S13). The minimization procedure is carried 
out numerically, with the main results for the magnetic order of the spin chain displayed in the exact phase diagram 
of Eig. 2 of the main text. 

We note that in the dilute limit ^ 1 we have GqT <C 1 and thus the off-diagonal components of are small. 
This is because Go rapidly falls off as thereby provides a controlled expansion for d > 1. We may then 

approximate the energy by expanding (S6) to second order in Go^- 

^=iE/ ^Tr[T,go,ijTjgoji]. 


(S15) 
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Carrying out the remaining trace and frequency integral leads to Eq. (6) of the main text in the limit r <C For 
r > ^ the terms in square brackets in Eq. (6) are modified by powers of r, but the entire expression oc becomes 

exponentially suppressed in this limit and we do not consider it further. As discussed in the main text, the perturbative 
expansion leading to the Eq. (S15) fails for < A/kpa due to strong YSR-Cooper pair hybridization. Formally, the 
l/l^l singularity in Eq. (6) is regularized by keeping all terms under the logarithm in Eq. (S6). In that case one must 
understand the structure of the term under the logarithm in Eq. (S13), which encodes the subgap YSR spectrum as 
well as its topological properties. 


YSR SUBGAP SPECTRUM 


The dispersion relation for the chain of YSR states is determined by having zero determinant under the logarithm 
in Eq. (S13). The sub-gap excitations induced by the spin lattice thus correspond to roots of the equation 


0 = det 


Qo{u,k - qa^/2) - T ^(w) 


(S16) 


where we have pulled out a factor of detT for convenience, since it has no sub-gap zeroes. One may also derive 
Eq. (S16) from the Dyson series for the exact electron Green function (see also Ref. [S3]), given by 


= + (S17) 

I 

where, within the spiral ansatz of Eq. (S8), we have Ji = ^ Multiplying Eq. (S26) on the left by 

^i^qaia^ and on the right by we obtain 


I 

In momentum space Eq. (S18) becomes, in the notation of the main text. 


(S18) 


Q{ikn, k - qcrz/2) = Go{ikn, k - ^ - q(Jz/2)JaxG{ikn, k - qcFzl2), (S19) 

while the inverse of Eq. (S19) may be written as 

G~^{ikn, k - qaz/2) = Go^{ikn, k - qaz/2) - Ja^. (S20) 

Splitting the bare Green function into non-local, and local parts 

Go{ikn,k) ^) + 0), (S21) 

j 


with Go defined in Eq. (S14), the local term and the self-energy may be combined into the T-matrix whereby 

G~^{ikn, k - qcFz/2) = Go{iK, k - qaz/2)~^ (T~^{ikn) - Go{ikn, k - qcfzl2^ Jcr^. (S22) 

After analytical continuation, ikn u i0+, the poles of G are found as zeroes of the determinant of G~^- The 

subgap solutions of Eq. (S21) occur solely from the zero determinant of the term in parentheses, leading to the same 

condition as Eq. (S16), which may be rewritten as 


0 = det 


B(7x + 


UJ + ATx 


{As AtCJz) ^CFz)'^z 5 


with B = {'kJv-f/2) ^ = ^^(5 + A)/(5 — A) and 


I = Rege{uj,k), a = Re^o(^,^), 

As = 1 + lmge{i^, k), At = Im^o(^, k). 


(S23) 


(S24) 
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Here the functions g^jo are defined as 


5e/o(w, k) = - 


1 

2kFCi 


ln(l — — ^-(a/C)y^l-uj^/A^+i(kF-k-\-q/2)a^j ^ ^ 

(S25) 


For a given k, Eq. (S23) has four roots placed symmetrically around zero, two of which lie in the continuum above 
the gap. The sub-gap excitation energy oj = generally takes a rather complicated form, which can be simplified 
by expanding to second order in and ^/A, to arrive at: 


Ek = ^(/ife-e)2+A2, (S26) 

where hk = Alm^e(0? k) and = —ARe^o (05 k). The identical spectrum Eq. (S26) was found in Ref. [S4]. We may 
now use the more general result (S23) to investigate the topological properties of the YSR chain. 


TOPOLOGICAL PHASE TRANSITION 

One may determine the topological classification by studying the exact inverse Green function at zero frequency, 
which defines a ‘topological Hamiltonian’ [S5] for the sub-gap electron degrees of freedom that propagate along the 
spin chain. As we show in Sec. , the precise form of the topological Hamiltonian is somewhat arbitrary and thus 
allows one to choose the most convenient representation. To that end, we define the topological Hamiltonian as the 
zero-frequency limit of the 4x4 matrix used in Eq. (S23) 

Htopik) = GoiO, k - - T“^(0) 

= Ba^; + Itz + aa^Tz + (S27) 

with coefficients defined as the zero-frequency limit of Eq. (S24). We have written i7top(^) in a suggestive form 
to make the analogy with the nanowire Hamiltonian of Refs. [S6, S7]. In this case, one would interpret B as an 
effective Zeeman field, A^/^ as singlet/triplet pairing, as kinetic energy and a as spin-orbit coupling. As discussed in 
Refs. [S8-S12], the topological class of i^top is BDI due to hidden time-reversal, and chiral symmetries, and is therefore 
characterized by a Z topological invariant. This can be shown by noting that, in addition to the usual particle-hole 
symmetry EHtop{k)E~^ = —Htop{—k) where E = (JyTylC and 1C denotes complex conjugation, there is an effective 
(anti-unitary) time-reversal operator O = —icFxlC such that 0HiQp{k)0~^ = i7top(~^)- The operator O satisfies 
0^ = 1 and is distinct from the canonical time-reversal operator 0 = icjylC^ 0^ = —1, which does not commute with 
i^top in the presence of the Zeeman field B. The existence of O allows one to define a chiral operator C = O 5 = (JzTy^ 
= 1, that anti-commutes with i^top- {C,i7top(^)} = 0- This implies that i7top(^) can be block off-diagonalized in 
the eigenbasis of C, a representation that may be achieved with the unitary transformation U = 

C/tfftopC/ = g ^ , (S28a) 

A = io. Ag + + (A^ -|- (S28b) 

The Z topological invariant measures the degree of the map from the Brillouin zone to the U(l) phase of the complex 
number z{k) = detA/|detA| , 7ri(5'^) = Z, and is therefore represented by a winding number W. Erom Eq. (S28b) we 
find 


{'kJpf/‘2)‘^ (1 - ijge - go)) (1 + ijgt + g^)) - 1 

\{nJpF/2r{l-i{g,-g,)){l + i{gt+g*))-l\ ^ ^ 

with 

w =E <i— = A i dkdk\nz(k). (S30) 

2t:i J z 27n J 

The value of W is quantized and can only change when 2 : becomes ill-defined. This occurs when detA = 0, and 
coincides with the closing of the spectral gap deti^top = 0, see Eq. (S28a). 
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Since the function go{^^k) vanishes at the time-reversal symmetric points k = 0,±7r/a, the location of the gap¬ 
closing transitions simplify in this case to give 


\l-igff-l 

where gf = ge{0, k) with A: = 0 or 7r/a, respectively. In the limit kpa ^ 1 we find 


(S31) 


— fs Imgf, 


(S32) 


which is equivalent to the gap closing transition at /c = 0, 7r/a in the YSR band for kpd ^ 1 (^/c = — s = Aj^ = 0)- 

Furthermore, in the limit ^ ^ a this leads to 

— ^ —(kpCi mod ^7r)//ci?a, (S33) 


see Fig. 2. Other gap closing transitions for k ^ O^n/a are not precluded and do occur. However, their locations in the 
phase diagram are difficult to determine analytically for arbitrary kpCi- An example of such behavior was discussed 
in the main text in the limit kpa ^ ^ and was shown to lead to the gapless ferromagnetic line (8). 


Topologically Equivalent Hamiltonians 


In this section we show that the definition of the Hamiltonian (S27) leads to the same topological classification as the 
Hamiltonian htop defined by the prescription of Ref. [S5] to be the negative inverse Green function at zero-frequency. 
According to Eq. (S22) this gives 


^top(^) — Qo{k 





J 1 


(S34) 


with the (zero) frequency variable suppressed. We may write htop = Hi^top^ = D — B (cf. Eq. (S20)), with D = Jgx 
and B = Go{k — Thus, if H, D anti-commute with the chiral operator C = (JzTy^ then Chio^{k)C~^ = —ht{k) 

has the same chiral symmetry as i^top- Since one readily sees CDC~^ = —D^ we need only prove that CBC~^ = —B. 
Since B~^ = i^top + clearly anti-commutes with C, one may simply take the inverse of the equation CB~^C~^ = 
—B~^ to prove that B also anti-commutes with C. Here we have utilized the fact that B is an invertible matrix for 
every k due to the bare electron spectrum being gapped by superconducting correlations (the zero eigenvalue may 
occur only in i^top)- 

In the eigenbasis of C, the matrices B also become purely off-diagonal, and the off-diagonal element of htop{k) is 
therefore of the form BAD^ where B and D are the off-diagonal components of B and respectively. The winding 
number W associated with phase of the complex number 


z'{k) 


detA. detH detC 
|detA| |detH| \detC\ 


(S35) 


is unchanged because the determinants of H, C never vanish, and hence W may only change when detA = 0. Thus, 
the topological classifications of htop and i^top are equivalent. 
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